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There is a growing interest in the relation between Bose-Einstein condensation (BEC) and su- 
perfluidity. A Bose system confined in random media such as porous glass is suitable for studying 
this relation because BEC and superfluidity can be suppressed and controlled in such a disordered 
environment. However, it is not clear how this relation is affected by disorder and there are few 
theoretical studies that can be quantitatively tested by experiment. In this work, we develop the 
dilute Bose gas model with a random potential that takes into account the pore size dependence 
of porous glass. Then we compare our model with the measured low-temperature specific heat, 
condensate density, and the superfluid density of 4 He in Vycor glass. This comparison uses no free 
parameters. We predict phenomena at low temperatures that have not yet been observed. First, the 
random potential causes a T-linear specific heat instead of the T 3 dependence that is usually caused 
by phonons. Second, the BEC can remain even when the superfluidity disappears at low densities. 
And third, the system makes a reentrant transition at low densities; that is, the superfluid phase 
changes to the normal phase again as the temperature is reduced. This reentrant transition is more 
likely to be observed when the strength of the random potential is increased. 



PACS numbers: 



I. INTRODUCTION 



Bose-Einstein condensation (BEC) and superfluidity 
of liquid 4 He in random environments including Aerogel 
and Vycor glass are active problems in quantum fluid re- 
search. In particular, finding out how spatial confinement 
affects the Bose fluid has stimulated both experimental 
and theoretical studies. 

Below the A temperature of 2.17 K, liquid 4 He enters 
the superfluid state and behaves as though it has no 
viscosity. Superfluidity is a macroscopic quantum phe- 
nomenon as well as superconductivity, and understand- 
ing both have been one of the major goals of quantum 
statistical physics. The various observations of superflu- 
idity was successfully explained by the phenomenological 
two fluid model 0] , which is based on the idea that the 
system consists of an inviscid superfluid and a viscous 
normal fluid. On the other hand, the lambda transi- 
tion had been thought to be caused by BEC, which was 
confirmed by neutron scattering experiments . With a 
BEC, a macroscopic number of particles occupies a single 
particle ground state and is described by a macroscopic 
wave function. The inviscid superflow can be described 
by this wave function g . However the relation between 
BEC and superfluidity is not completely understood. Al- 
though superfluidity and BEC are closely related to each 
other, one is not necessary or sufficient for the other. For 
example, in a two-dimensional Bose system, Kosterlitz 
and Thouless proved that superfluidity can exist even 
without BEC H , and the superfluidity was actually ob- 
served in 4 He films || . The Bose system in a random 
environment might be another good example for study- 
ing the relation between BEC and superfluidity jfj, Q . 
This system has received considerable attention because 
localization effects allow condensed particles to belong 
to the normal fluid rather than the superfluid and BEC 



can separate from the superfluid. The phase diagram of 
this system has been discussed, showing a specific non- 
superfluid phase. Thus, studies of this system can reveal 
the relation between BEC and superfluidity. 

Porous glass such as Vycor is often used as a random 
media in experimental studies. Vycor glass is a from 
30% to 70% porous glass containing wormholelike pores, 
the characteristic diameters of which vary from 30A to 
100A. By adjusting the pore size and the adsorbed 4 He 
coverage, we can change the density of He and the super- 
fluid transition. By using torsional oscillators, Reppy et 
al. H, ^ interesting features at various pore sizes of Vycor 
glass or coverages were observed, particularly the behav- 
ior of the superfluid critical temperature and the tem- 
perature dependence of the superfluid density. The su- 
perfluid component has two-dimensional behavior when 
the pore size of Vycor glass is large, and becomes three- 
dimensional as the pore size is reduced. The superfluid 
density in such porous glass is smaller than that of bulk 
4 He and its critical temperature decreases with the cov- 
erage. Below a certain coverage, the superfluid density 
can no longer exist, even near K. These results show 
that superfluidity is broken by the random environment. 
It is also important to find out how disorder affects the 
BEC. BEC and its elementary excitation in liquid 4 He 
can be observed by neutron scattering. Bulk 4 He has the 
usual excitations such as phonons, maxons and rotons [jlj. 
They are collective excitations on BEC, the excitations 
except for phonons being absent above the critical tem- 
perature. Dimeo et al. jlCj] and Plantevin et al. [[fl) used 
neutron scattering and a torsional oscillator to measure 
the elementary excitations and the superfluid transition, 
respectively, of 4 He in porous glass. Surprisingly, the dis- 
persion curve in porous glass was the same as that in the 
bulk, which means that the disorder does not affect the 
elementary excitations. Furthermore, these elementary 
excitations were observed even above the superfluid crit- 
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ical temperature. Hence, the BEC might persist in the 
disordered system even above the superfluid critical tem- 
perature. These results have not been understood com- 
pletely yet; thus, the exact relationship between BEC 
and superfluidity remains puzzling. 

This problem is also interesting from the following the- 
oretical standpoint. In a Bose system confined in a ran- 
dom media, the long-range-order correlation due to BEC 
can compete with the disorder, so that the BEC critical 
temperature can be reduced. Huang and Meng proposed 
a model of the three-dimensional dilute Bose gas in a 
random potential |L3) that assumed a small coverage of 
4 He in Vycor glass . Because it is difficult to formu- 
late the random potential for the porous glass, they used 
a delta-functional impurity potential and analyzed their 
model using the Bogoliubov transformation and taking 
an ensemble average. They found that both BEC and 
superfluidity are depressed by the random potential and 
the superfluidity disappears below a critical density, even 
at K, which is qualitatively consistent with the obser- 
vations by Reppy et al. They also predicted a reentrant 
transition at low densities; that is, the superfluid phase 
enters a normal phase again with decreasing tempera- 
tures. However, the random potential of their model does 
not include the pore size, and thus it is difficult to quan- 
titatively compare to experimental results for a range of 
pore sizes. Another model is the Bose Hubbard model 
with the random potential. By considering the transfer 
energy, the on-site repulsion, and the random potential, 
Fisher et al. jl4| found that the Bose glass phase can 
exist with the superfluid phase and the Mott insulating 
phase. The Bose glass phase is similar to the Anderson 
insulating phase [051 in metal. In the Bose glass phase, 
the condensed particles are localized and thus do not con- 
tribute to superfluidity. Thus, the Bose glass phase could 
influence the collective excitations even above the super- 
fluid critical t emp erature. However, the theoretical exci- 
tation energy for the Bose glass phase disagrees with 
measurements by neutron scattering experiment Jll| , so 
it is not yet clear whether the Bose glass phase has ac- 
tually been detected. Finally, it should be noted that 
Huang and Meng's model cannot describe the Bose glass 
phase because the ensemble average makes the system 
uniform. 

Few theoretical studies of this random system are 
quantitatively compare to the experiment. Thus, in 
this work, we improve Huang and Meng's model |l3[ | by 
adding the size dependence of the random potential in- 
stead of using their the delta- functional potentials. The 
strength of the random potential can be estimated by 
comparing calculated and experimental critical coverages 
below which the superfluid density disappears, even at 
K. As a result, our model has no free parameters and 
can be used for quantitative comparisons to experimen- 
tal data. This enables us to determine whether or not 
our picture of the three-dimensional dilute Bose gas in 
random potential is applicable to a real system. Our for- 
mulation cannot address this question at high tempera- 



tures due to the high number of thermally excited quasi 
particles. As far as the condensate density is almost in- 
dependent of temperature at low temperatures, however, 
our formulation works well, leading to the following re- 
sults. (1) The specific heat agrees quantitatively with ex- 
perimental data at low temperatures. (2) Because of the 
random potential, the specific heat is not proportional to 
as occurs for phonons, but to T . Furthermore, by 
obtaining the condensate density and the superfluid den- 
sity, we found the following. (3) When the total density 
is sufficiently low, BEC can persist even when the su- 
perfluid density disappears below that critical coverage. 
(4) The random potential causes a reentrant transition 
of the superfluid phase. Finally, we show why decreasing 
the open pore density of the Vycor glass should allow the 
reentrant phase to be detected experimentally 

A brief summary of our paper is as follows. In Sec. II, 
we describe our model of the dilute Bose gas in a random 
potential and derive the partition function. Section III 
tests our model by quantitatively comparing calculated 
to experimental specific heats. In Sec. IV, the BEC 
density and the superfluid density are obtained and their 
characteristics are discussed. Section V is discussion and 
conclusions. 



II. MODEL 

Superfluid 4 He adsorbed in Vycor glass can be mod- 
eled by a three-dimensional dilute Bose gas in a random 
external potential |13| . The grand canonical Hamiltonian 



H - fiN = K 



— V 2 + U{x)-n V(x) 
2m 



+ 



f'o 



/rf. W i'HiHi(4 (i) 



where ^(x) is the field operator for Bose particles of 
mass m, N = J dx 3 ^^ (x)'$>(x) is the number opera- 
tor, /i is the chemical potential, and U(x) is the exter- 
nal random potential that represents the effect of Vycor 
glass. The first term of the Hamiltonian is the kinetic 
energy and the external potential, whereas the second 
term refers to the hard-sphere interaction between par- 
ticles with vq — Ait ah 2 /m being the coupling constant 
with the s-wave scattering length a. This repulsive in- 
teraction prevents all particles from being localized at 
the minimum of U(x). This has similarities to the Fermi 
system with disorder [Tt], |l8); for example, Fermions 
cannot localize in a single orbital in space due to the 
Pauli exclusion principle. Therefore, the Fermion system 
is stable even if it is free from the repulsive interaction. 
On the other hand, to prevent the system from collaps- 
ing into the minimum of U(x), the Bose system should 
include a repulsive interaction. This makes the problem 
more complicated than that of the Fermi system. 
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Proceeding in a standard fashion, we introduce the 
free particle annihilation and creation operators and 

&L- We assume that the level with k = is macro- 
scopically occupied with occupation number Nq, so clq 



and a are replaced by a c-number y/No. By making 
a Fourier transformation and neglecting all off-diagonal 
terms Ut,a\,au'' and VoaXa\,au''auiii, we obtain 
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where V is the volume of the system, no = Nq/V is 
the number density of condensate, and Uy, is the Fourier 
transformation of U(x). By neglecting the off diagonal 
terms, we are neglecting the interactions between the ex- 
cited particle and the random potential and that between 
pairs of excited particles; these become important as the 
temperature rises and the condensate density decreases. 
Hence, this approximation is poor when many particles 
are thermally excited. Nevertheless, this approximation 
is useful at low temperatures where the condensate den- 
sity is almost independent of temperature. All results 
here are obtained for these low temperatures. 

This Hamiltonian can be diagonalized by the Bogoli- 
ubov transformation 
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Then the coefficients 7^, and the quasiparticle spec- 
trum ujfe are given by 
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(4b) 
(4c) 
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Next, we take an ensemble-average to quench the ran- 
dom potential. The random potential simulates Vycor 
glass with a characteristic pore size r p as follows. The 
quenched potential Uu may decay above the character- 
istic wave number k p — 2n/r p . Thus, we assume the 
averaged potential 
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( U k U -k)^ = R o ex P 



2k 2 
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where av denotes the ensemble-average. i?o, with dimen- 
sion (energy) 2 (length) 3 , is the characteristic strength of 
the random potential. Equation ([5]) makes our model 



completely different from Huang and Meng's, and we will 
show that the results are also different. The coherence 
length of the BEC is thought to be from hundreds to 
thousands of A, whereas the spatial scale of disorder is 
the pore size in the glass, which is dozens of A. Hence, 
the macroscopic wave function of BEC is not sensitive 
to disorder in and between pores but instead depend on 
the disorder averaged over the coherence length. Hence, 
the ensemble-averaged system can become nearly uni- 
form. For a uniform Bose system, it has been proven that 
the elementary excitation spectrum becomes the gapless 
Goldstone mode fl9fl . Thus, we set A = in Eq. (Ed). 

The resultant diagonalized and ensemble-averaged 
Hamiltonian is 
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(6a) 

(6b) 
(6c) 

(6d) 
(6e) 



where t\ is the hard sphere interaction energy at K, 
similarly, 6r is that for the random potential. The quasi- 
particle spectrum u>k is the same as that in the hard 
sphere Bose gas model and is independent of the ran- 
dom potential. This independence is confirmed by neu- 
tron scattering experiments which justifies the above 
assumption of A = 0; conversely, if A ^ 0, the spectrum 
would depend on the random potential. 

This Hamiltonian enables us to obtain the grand par- 
tition function Q = Tr{exp(— /3K)} and various physical 
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quantities. The condensate density is defined by the fol- 
lowing relation: 



no = n 



fc^o 



where n is the particle number density. The second term 
represents the noncondensate particle number as 
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Here ni is the noncondensate density excited by the hard 
sphere interaction, tir is the density due to the scattering 
of condensate particles with the random potential, and A 
is the thermal de Broglie wave length. When a vanishes, 
tir becomes infinite. This means that the system would 
collapse if there were no repulsive interactions between 
particles. 

Because superfluidity is described by the two fluid 
model, the particle density n consists of the normal fluid 
density n n and the superfluid density n s . The super- 
fluid density n s can be calculated by linear response 
theory |20| . Because of its viscosity, only the normal fluid 
responds to a small, applied velocity field. Thus the nor- 
mal fluid density can be defined by the response of the 
momentum density ji(x,t) to the external velocity field 
Vi(x,t). Linear response theory gives the following rela- 
tion: 



ji(x,t) = Xij(x,t)vj(x,t), 
Xij(x,t) = ([j 4 (:z,i), .7,(0,0)]}, 
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where ^>(x,t) is the Heisenberg field operator. The static 
susceptibility Xy(fe) is defined as 



du d 3 k ,t ;fc t . 
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Because of the rotational invariance, the static suscepti- 
bility Xij (k) can be written 



k>i kj 



B(k), 
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where A(k) and B{k) are the longitudinal and transverse 
parts, respectively. The transverse susceptibility B(0) is 
the normal fluid mass density. The superfluid number 
density n s is n — B(Q)/m. The susceptibility B(0) can be 
calculated by the Bogoliubov transformation in Eq. (||). 
After some tedious calculations, the resultant superfluid 
density is given by 



n s =n - n n i - n nR , 



dt- 



(1-e- 



n n R = ~n R , 



(12a) 
(12b) 

(12c) 



where n n \ is the normal fluid density due to the ele- 
mentary excitations, and n nR is that due to scattering 
with the random potential. The density n n \ can be also 
obtained using Khalatnikov's method that is based on 
Galilean invariance pl| . The relation n n R — A/3ur — 
71 r + l/3riR shows that the random potential causes the 
larger normal fluid density than the noncondensate den- 
sity; some condensate particles are captured by the ran- 
dom potential to participate in the normal fluid. This 
makes it possible to destroy superfluidity even at K 
when n n R becomes comparable to n. 

This formulation can be used to obtain various phys- 
ical quantities including the condensate density, the su- 
perfluid density, and the specific heat. 



III. COMPARISON WITH EXPERIMENTS 

In this section, we compare the calculated specific 
heat and the superfluid density with experimental re- 
sults. Quantitative agreement is shown to be good at 
low temperatures, which supports our assumption of a 
dilute Bose gas in the random potential. Furthermore, 
we show that the random potential leads to as-yet unob- 
served behavior of the specific heat. 

To make a quantitative comparison, we give the fol- 
lowing numerical values to the parameters: m ~ 6.6 x 
10~ 27 kg and a ~ 5x 10~ 10 m are the mass and the s-wave 
scattering length of a 4 He atom. Other parameters are 
from the experiments of Reppy et al. [||. The volume V 
of open pores in the Vycor glass (about 40% of the total 
volume of the Vycor glass) is about 1 cm 3 . The parti- 
cle density n of 4 He inside the Vycor glass is estimated 
as follows. In Vycor glass, the atoms are adsorbed and 
fully cover the surfaces of the open pores due to the van 
der waals attraction. The pore area is about 10 8 m 2 /m 3 . 
The rest of atoms, which do not participate in the first- 
layer solid, can behave as a dilute gas inside the pores. 
The particle density n of the dilute gas is obtained by 
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subtracting the adsorbed amount from the total amount. 
This density is estimated to be from 0.001% to 70% of 
the density of bulk liquid 4 He Ubuik ~ 2.1 x 10 28 /m 3 . 
Because the first layer of 4 He adsorbed on the surfaces 
cannot move and behaves as a solid, we assume the pore 
size 30A of Vycor glass is effectively reduced by 2a. Thus, 
r p is estimated to be 20 A. The last parameter Po, which 
is the strength of the random potential, can be fixed by 
comparing to experiment. Shown in Fig. |l|(a) are the 
data of zero temperature superffuid signals taken in an 
experiment that used a torsional oscillator (Fig. 12 of 
reference Q). Because the superfluid component docs 
not contribute to the moment of inertia, the resonant 
frequency and the period of oscillation differ from those 
without superfluid. The period difference AP is approx- 
imately proportional to the superfluid component. Here, 
the superfluid density is nearly proportional to AP and 
disappears at a coverage of 17.5 mg. Figure |l|(b) shows 
the superfluid density at K from Eq. (|l2]). As in the 
experiment, the superfluid density becomes zero at a cer- 
tain coverage that depends on P . Thus, the value of 
Po can be fixed using the comparison with Fi g. p]( a); i.e., 
Ro = 5x 10 -75 J 2 m 3 . Here, we define R w = y/Ron, which 
is the single particle energy converted from Po. In the 
Vycor glass, R w /kB is about 0.001 ~ IK. Just above the 
critical coverage, the superfluid density increases linearly 
for both the experiment and the calculation; however, 
their slopes cannot be compared because the amplitude 
of n s is unknown in the experiment. 

Because all parameters are now fixed, we will quantita- 
tively compare calculations to experiments. The specific 
heat can be obtained from temperature differentiation of 
the free energy: 



fir 
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where fii is the free energy of the elementary excitation 
and the hard sphere interaction, and fi/j is the free en- 
ergy from the random potential. Figure |^ compares our 
results to data of low temperature specific heat taken 
from Fig. 1 of reference In Fig. ||(a), which shows the 
data for high density, the density n is fixed from the ex- 
perimental coverage, whereas we fix the density from the 
superfluid critical temperature in Fig. |^(b) at low den- 
sity. This is because we have no information about the 
data on the coverage. Figure |^(b) also shows the super- 
fluid density. The theoretical results agree quantitatively 
with experiment without using free parameters. Above 
1.0 K in Fig. ||(a), the calculated condensate density be- 
gins to decrease rapidly; here, our criterion of constant 
condensate density fails, which likely causes the discrep- 
ancy with experiment. However, Fig. |^(b) shows that 
the calculated specific heat agrees with experiment up to 
temperatures near the superfluid critical temperature; in 
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FIG. 1: Superfluid signals of experiments (a) and calcula- 
tions (b) near K. AP is the resonant period difference in a 
torsional balance experiment, which is approximately propor- 
tional to the superfluid density M . 



this temperature region, the calculated condensate den- 
sity hardly decreases. This means that the system is 
more dilute than that of Fig. ||(a) and thus is affected by 
the random potential rather than the elementary excita- 
tions. We discuss this in the next section. Nevertheless, 
these comparisons show that our model is accurate at low 
temperatures. 

Our model predicts an effect at low temperature that 
is due to the random potential but has not yet been ob- 
served. This is shown in Fig. ||, which is the log-log plot 
of Fig. ||(a). When the system is free from the random 
potential, the specific heat should increase a T 3 because 
of the contribution from the phonons. However, with a 
random potential, the dependence is linear in T at low 
temperatures. This means that at low temperatures, the 
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FIG. 2: The specific heat data from in experiments (plot) 
and calculations (line). In (a), experimental data are given 
by Fig.l of reference ||. The circles, triangles and squares 
respectively correspond to full pores (a = 1), a = 0.780 and 
a — 0.636. Here a is the ratio of the coverage to the full 
pore coverage. In (b), calculated and experimental superfluid 
densities are compared (T c = 0.163K). 



contribution from the random potential is larger than 
that from the elementary excitations (phonons) . The free 
energy Qr in Eq. (|l^) depends on the temperature only 
through the condensate density uq. Qr is the energy 
from the scattering of the condensate particles with the 
random potential, and the resultant specific heat is given 
by the energy that the condensate particles need to slip 
out of the random potential. Experimental observation 
of this T-linear dependence might clearly identify the in- 
fluence from the random potential. 



FIG. 3: The log-log plot of Fig. |(a). Data of the specific 
heat at Ro = 0(n/ni,uik = 0.35). Two lines for tx T and cx T 3 
are added. 



IV. CONDENSATE DENSITY AND 
SUPERFLUID DENSITY 

This section describes some characteristic behavior of 
the condensate density no and the superfluid density n s 
derived from our model. 

Figure. |J shows the dependence of no and n s on tem- 
perature and density. Both no and n s decrease with de- 
creasing density, even at K. This means that the effect 
of the random potential on no and n s becomes larger as 
the density is reduced. Figure. |] shows clearly the differ- 
ence between n s and no- Below the critical density, the 
superfluid density disappears, although the condensate 
persists. This situation indicates that the condensate 
particles cannot move as a superfluid because they are 
trapped by the random potential. We expect that this 
theoretical result will be confirmed by measurements of 
condensate density. 

Figure. || shows the temperature dependence of the 
superfluid density just before superfluidity disappears. 
This figure shows the reentrant transition at which the 
superfluid density n s goes to zero with a decrease of tem- 
perature. In this temperature region, the condensate 
density is almost constant; hence, our formulation should 
work well in accordance with the criterion described in 
Sec. II. This reentrant transition is understood as follows. 
The condensate depletion n r of Eq. (||c) comes from the 
scattering of the condensate particles with the random 
potential. This decreases with hq as the temperature 
rises, so that n n R = 4/3n^, of Eq. (p~2]c) also decreases. 
As shown in Fig. ^|. the magnitude of this decrease in 
n n R exceeds the increase in n„i of Eq. (JT^b) , which is 
the normal fluid density due to the elementary excita- 
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(a) 




(b) 



FIG. 4: Temperature and density dependence of rio (a) and 
n a (b) at R = 5.0 x l(T 75 J 2 m 3 . 



tions, in the very low temperature region T < 0.5 mK. 
In other words, condensate particles that are trapped by 
the random potential at lower temperatures can escape 
at higher temperatures and thus participate in superflu- 
idity. This reentrant transition has not been observed 
experimentally, probably because it should only occur at 
very low densities and low temperatures. However, large 
values of i?o can make the reentrant transition observ- 
able as follows. We define the temperature T max as that 
which maximizes the superfluid density n s , and define 
An s = n s (T = T max ) - n s {T = 0). Fig. | shows the 
density dependence of T max and An s . Both variables 
increase with Rq. Therefore, the reentrant transition is 
more likely to be measured at larger Rq. The parameter 
Rq is the strength of the random potential over the entire 
space, and one way to increase Rq is to decrease the open 
pore density of the Vycor glass. 




"6.30 x 10 n 



"6.25 x 10 n 



0.5 1 1.5 

Temperature[mK] 



FIG. 5: Temperature dependence of n s at low temperature 
and low density. In this regime, the superfluid density n s goes 
to zero with a decrease of temperature (reentrant transition). 




0.5 1 1.5 2 2.5 3 3.5 4 
Temperature[mK] 

FIG. 6: Temperature dependence of n„i and n n R near the 
reentrant distribution at n — 6.25 x 10~~ J ribuik (dashed line 
in Fig. R). 



V. CONCLUSIONS 

The present paper describes the dilute Bose gas system 
in a random potential. The outcomes of our studies are 
as follows. 

By including the pore size dependence of Vycor glass in 
the random potential, our model could closely match the 
experimental conditions of liquid 4 He in Vycor glass. We 
fixed the strength of the random potential by equating 
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the theoretical and experimental critical coverages below 
which the supcrfluid density at K vanishes. No other 
parameters could be adjusted, and thus we could quanti- 
tatively compare theory to experiment for other physical 
quantities. 

First, we showed that the calculated specific heat for 
Vycor glass quantitatively agrees with measurements. 
This indicates that liquid 4 He in Vycor glass behaves as a 
dilute Bosc gas in a random potential. For low tempera- 
tures, the calculated specific heat was linear in T because 
of the random potential. Second, the BEC was shown 
to persist even when superfluidity disappears below the 
critical density. Finally, we showed that the reentrant 
transition of the superfluid phase is more likely to be ob- 
served experimentally by increasing the strength of the 
random potential. 

Because we neglected interactions between pairs of ex- 
cited particles and between excited particles and the ran- 
dom potential, this model does not apply to systems at 
high temperatures. To overcome this limitation, we are 
improving the model to include these interactions and 
will report on this more general model in the near fu- 
ture. 
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